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A phenomena-logical analysis of the scalar meson /o(980) is performed that relies on the 
quasi- two body decays D and D s — > /o(980)P, with P = tt, K. The two-body branching 
ratios are deduced from experimental data on D or D s — > nmr, KKn and from the 
/o(980) — > 7r + 7r~ and /q(980) — » K + K~ branching fractions. Within a covariant quark 
model, the scalar form factors for the transitions D and D s — > /o(980) are computed. The 
weak D decay amplitudes, in which these form factors enter, are obtained in the naive 
factorization approach assuming a qq state for the scalar and pseudoscalar mesons. They 
allow to extract information on the Jo (980) wave function in terms of uu, dd and ss pairs 
as well as on the mixing angle between the strange and nonstrange components. The weak 
transition form factors are modeled by the one-loop triangular diagram using two different 
relativistic approaches: covariant light-front dynamics and dispersion relations. We use the 
information found on the /o(980) structure to evaluate the scalar and vector form factors in 
the transitions D and D s — > /q(980), as well as to make predictions for B and B s — > /q(980), 
for the entire kinematicahy allowed momentum range of q 2 . 



I. INTRODUCTION 



Scalar mesons have been a recurrent topic over the past 30-40 years. Whereas the existence 
of the cr//o(600) has been a longstanding open question since the 1960s, the /o(980) and ao(980) 
were firmly established in tttt scattering experiments in the 1970's The known ++ mesons fall 
into two classes: near and about 1 GeV and in the region 1.3—1.5 GeV. The scalar objects below 
1 GeV form an SU(3) flavor nonet [2(. This nonet contains two isosinglets, an isotriplet and two 
strange isodoublets. Among these lighter scalars, the isosinglet /q(980) and the isotriplet ao(980) 
are rather narrow with their widths T ranging from 40 to 100 MeV [3j. Both scalars strongly 
couple to the KK channel and lie close to its threshold at 987 MeV. This closeness alters the 



shape of their resonant structure and the description of the /o and ao requires a coupled-channel 



scattering analysis. The simple quark model views these scalar mesons as orbitally {L = 1) 
excited qq states and has been advocated, for example, by Tornqvist and Roos [4] as well as in 



Ref. fl. 



However, some studies 



tend to favor four-quark configurations of the scalar mesons, 



as do coupled-channel analyses |7j or potential models of molecular states strongly coupled to 
7T7T and KK channels [g]. 



[9 



The emergence of the /o(980) as a pole of the tttt amplitude in the S wave [9J is also 
well established in three-body decays of B mesons Recent tttt effective mass range 

distributions, obtained from an isobar model fit of the B — ► 7r + 7r~ K and B —* KKK Dalitz 



plots by the Belle 



n 



1.3] and BaBar Collaborations Q, Q, [3], display distinct peaks 
about 1 GeV. Scalar resonances have also been observed in the charmed three-body decays 



D -► TrTTTr-TTTrK, KKK at CLEO 



E687 



24, 



0KFOCUS Q 



251 ] . Remarkably, in Ref. 



FOCUS 1181, ARGUS [191, BaBar [201, E791 



21 



22 



[231 ] an experimental evidence for a light and broad scalar 
resonance in the m nn spectrum of the D — > tttttt decay was found, which may be identified with 
the /o(600) and a peak within the /o(980) mass range is also observed. Although a considerable 
amount of data has been accumulated over the years, it has yet not been possible to elucidate 
the precise _/o(980) quark structure, i.e. whether one deals with a two-quark or rather a four- 
quark composite, and thus far there is no consensus on that matter. On the other hand, viewing 
the /o(980) exclusively as a qq or q 2 q 2 state may simply be too naive 26|]. In this context, an 



interesting proposition to shed light on the constituent composition of the /o(980) was recently 
made by Maiani, Polosa and Riquer [271 ] . Their method consists in comparing the ratio of the 
decay rates Df — ► tt + (K + K~) and — ► it + {KsKs). This ratio is predicted to be 1/2 if the 
/o(980) is an / = 0, qq state, whereas the composition Jq = [sq][sq],q = u,d, could yield a 
different value owing to possible interference patterns between 1 = and 1 = 1 amplitudes in 
the tetraquark picture of these decays. For a general overview on scalar mesons, we refer to the 
Particle Data Group review Q] and references therein. 

In the case of B — > fo(980)K decays, one may advance plausible reasons to limit oneself to the 
qq picture of the /o(980). Because of the large B mass, the outgoing mesons are virtually massless 
particles, which prompts to expand the corresponding bound states in terms of Fock states. 
Quark configurations like q 2 q 2 or q 2 q 2 q therefore belong to higher Fock states. A handwaving 



argument by Cheng, Chua and Yang 



281 ] suggests that the qq component of the energetic /o(980) 



may be more important, as two rapid qq pairs are less likely to form a fast moving /o(980). In 
our models we neglect higher Fock contributions to the /o(980) bound state. 



In the two-quark model of the light scalar octet below 1 GeV, assuming an ideal SU(3) mixing 
angle of the octet states, the flavor content of the fo (980) is purely strange (ss) while that of 
the a or /o(600) is purely nonstrange [uu + dd) (see e.g. [28]). In such a picture the a is the 
lightest scalar, the /o(980) the heaviest. However, there is compelling experimental evidence, 
for instance from _/b(980) — » nir decays Q], that the /o(980) cannot be made of strange quarks 
only. We therefore introduce in this work some mixing between the strange and nonstrange 
flavor content. Experimental implications on this mixing have been the object of several studies 



(see e.g. 



29, 



2, 



31 



22l) 



In this paper we complete preliminary calculations of B — > /o(980) and D —* /o(980) 
pseudoscalar to scalar (P — > S) transition form factors [33I [3^ . Transition form factors are 
important for an understanding of the hadronic component of heavy-to-light decay amplitudes 
and their precise evaluation is crucial to a reliable determination of Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements. The short-distance physics is calculated in perturbative QCD, which 
comprises radiative vertex corrections to local four-quark operators in the operator product 
expansion 35f] as well as hard-scattering corrections with the spectator quark that go beyond 
the leading order 3f|- In contrast, the transition form factors are by nature long-distance 
nonperturbative hadronic matrix elements. They provide one ingredient of the factorizable 
amplitudes of the nonleptonic B decays mentioned above. A variety of theoretical approaches 
to heavy-to-light transition (pseudoscalar to pseudoscalar) amplitudes exist, either using light- 



cone sum rules [37J , light-front 



39 



40 



41 



42( | . Most recently, a 



or relativistic quark models 
comprehensive set of B meson heavy-to-light transition form factors, calculated with truncated 
amplitudes based on Dyson-Schwinger equations in QCD, was reported in Ref. (4^. Whereas 
the methods of Refs. 



37 



transfers q 2 , those in Refs. 
momenta. To our know 
with QCD sum rules 



381 only 



39, 



40, 



provide form factors for a small domain of timelike momentum 



41 



,142, 



28, 



431 ] give access to the entire range of physical timelike 
edge, P — > S transition form factors have only been evaluated so far 
441 ] at q 2 = four-momentum transfer. A functional extrapolation is 



required to access all timelike q 2 in these studies. 

The present work relies on two explicitly covariant formalisms: the covariant light-front 
dynamics (CLFD) and the dispersion relation (DR) approaches. Both require two size 
parameters as well as a mixing angle between the uu, dd and ss components to specify the 
/o(980) wave function. In order to deduce these parameters from experiment, we fit D and 
D s — * /o(980)P branching fractions, where P can be a pion or a kaon. Initially, decay amplitudes 
at tree level in the naive factorization approach are employed and neither annihilation nor 



4 

penguin topologies are considered. Already at tree level, nonleptonic two-body D decays can 
be reasonably reproduced within this simple factorization since penguin amplitudes are strongly 
CKM suppressed. However, since the charm mass m c is lighter than the bottom mass by roughly 
a factor three, nonperturbative contributions of order Aqcd/^c are more important than in B 
decay amplitudes. The factorization approach may then be less reliable. In order to study 
the discrepancy between theoretical and experimental branching fractions, we study the effect 
of phenomenological annihilation as well as penguin amplitudes. The decay amplitudes are 
proportional to the D and D s — ► /o(980) transition form factors we are interested in. 

The paper is organized as follows. In Sec. HU we present the CLFD formalism and give a brief 
review on the DR approach. The scalar /o(980) bound-state structure is described in Sec. lIIII In 
Sec. IIVt we list all physical constraints imposed in our model, namely experimental D-branching 
ratios and wave-function normalizations. The electroweak decay amplitudes, the D-decay tree 
topologies and all numerical inputs needed are presented in Sec.[V] In Sec. I VII we introduce the 
P — > S transition form factors, derived in CLFD and DR approaches. Details about the initial 
pseudoscalar wave functions, the pseudoscalar decay constant in the constituent quark model 
and the calculation of the P — > S transition form factors are given in Appendices [A] and [Bl 
In Sec. IVII| we discuss the fitting method, give numerical results for the theoretical branching 
ratios and then compare — ► /o(980) and Br s -\ — > /o(980) transition form factors obtained 
in both relativistic approaches 1 . The final Sec. IVIIll summarizes our work and some conclusions 
are drawn. 

II. TWO DIFFERENT RELATIVISTIC FORMALISMS 
A. Covariant Light-Front Dynamics 

In CLFD [451 ], the state vector which describes the physical bound state is defined on the 
light-front plane given by w-r = a, where u is an unspecified lightlike four vector (uj 2 = 0) 
which defines the position of the light-front plane and r is a four vector position of the system. 
CLFD proposes a formulation in which the evolution for a given system is expressed in terms 
of covariant expressions. Any four vector describing a phenomenon can be transformed from 
a system of reference to another by using a unique standard matrix which depends only on 
kinematical parameters and on u. The particle is described by a wave function expressed in 



1 Here and in the following, the notation refers either to the D or to the D s meson and similarly Br 3 -\ refers 
to the B or to the B B meson. 



terms of Fock components of the state vector which respects the properties required under any 
transformation. The meson of mass M will be described as a bound state of two constituent 
quarks with four momenta k\ and ki . The state vector describing this meson of four- momentum 
p, defined on a light-front plane characterized by u, is given by: 

\p,X) u = (2vr) 3 / 2 | *J A ffliaOJ| (fc 1) te jJ ,, W r)ot 1 (ki)at iI (k 2 )|0) 

x jWfo + k 2 -p - urr) exp(^)2(^)dT {2 J/^ (D 



where = y k? + m? and kj is the momentum of the quark i with mass mj . The parameter r 
is entirely determined by the on-mass shell condition for the individual constituents. In Eq. dl|) 
A is the projection of the total angular momentum J of the system on the z-axis in the rest frame 
and <7j is the spin projection of the quark i in the corresponding rest systems. We emphasize 
that the bound state wave function is always an off-energy shell object, with due to the 

binding energy, and depends on the light-front orientation. From the delta function ensuring 
momentum conservation, one gets: 

V =p + ujT = k 1 + k 2 . (2) 

To keep track of this conservation law, a momentum, ujt, is assigned to the spurion but there 
is no fictitious particle in the physical state vector, (see Fig. 1). The two-body wave function 
k 2 ,p, wt) in Eq. ([1]) can be parametrized in terms of various sets of variables. In order 
to make a close connection to the nonrelativistic case, it is more convenient to introduce the 
following pair of variables 45[ defined by 



k^-iCP^-^^-^-j, (3) 

\L i (V)u\ LO • p 

where V = k\ + k 2 , and L^ 1 ^) is the Lorentz boost. The momentum, k, corresponds, in 
the center of mass frame where ki + k2 = 0, to the usual relative momentum between the 
two particles. Note that this choice of variable does not assume that one is restricted to this 
particular frame. The unit vector n corresponds, in this frame, to the spatial direction of a;. 

One introduces the variables x and the vector R\ = (Rq, Rj_, Rii) where Rj_, Rn denotes the 
perpendicular and parallel components to the direction of the light-front: 

x = -, R\ = k\ — xp. (5) 

L0-p 



Since by construction Ry-uj = 0, and thus R 2 = — R^_, the light-front coordinates, which one 
will use in the present work, are then (x, Rj_). These variables can be expressed in terms of the 
ones in Eqs. ([3]) and dH). All details can be found in Ref. [451 ] . 

In terms of the variables (x, Rj_), we have for the relative momentum between two quarks of 
different masses: 

|r^_ + [(x — l)mi — X7712] 2 j |r^_ + [(x — l)mi — X7712] 2 j 



k 2 



(6) 



4x(l — x) R^ + (1 — x)m\ + xm\ 
B. Dispersion Relation approach 

The dispersion relation approach, in the context of the relativistic quark model, leads to 
transition amplitudes expressed as relativistic spectral integrals over spectral densities of the 
corresponding Feynman diagrams. Here we closely follow the derivation of Melikhov 42] 



to calculate the P — > S transition form factors. These are given by the double spectral 
representation over the square of the invariant masses of the initial and final quark-antiquark 
bound states. The spectral functions involve the wave functions of the participating mesons and 
the double discontinuities of the corresponding triangle Feynman diagram. Use of the Landau- 
Cutkosky rules allows to calculate these discontinuities and hence the transition form factors 
in the space- like region q 2 < 0. An analytical continuation in q 2 gives the form factors in the 
timelike region q 2 > 0. 

As in Sec. Ill A( the meson of mass M is a bound state of two constituent quarks of mass m\ 
and m,2 and four-momentum k\ and ki with 

s = (k\ + k^f \ k\ = m 2 , k\ = m|. (7) 

The relativistic bound state corresponds to a pole in the amplitude at s = M 2 and one can 
define a bound state wave function ip(s) in the vicinity of the pole by 

The function G v (s) in Eq. ([5]) represents the vertex of the bound state transition to the 
constituent quarks. The constituent-quark rescatterings lead to the normalization condition [2] 

f° p{s,m 1 ,m 2 ) ^ , . 

y (mi+m2)2 <s ~ m 2 ) 2 - i; w 

where the spectral density p(s,mi,m2) for a pseudoscalar meson reads 



pp(s, m 1 ,m 2 ) = — [s - (mi - m 2 ) J V{s - {mi + m 2 ) J , (10) 



7 

while for a scalar meson one has 

ps{s, m 1 ,m 2 ) = [s - (mi + m 2 ) 6{s - (mi + m 2 ) J . (11) 

87TS 7 

In Eqs. (fTUj) and (fTT|) . A(s,m 2 ,m2) is defined as 

A(s, m 2 , m|) = (s + m\ — m, 2 .) 2 — 4sm 2 , (12) 

and 6{z) is the step function, 9{z) = 1 for z > and 9{z) = for z < 0. 

Prom Eqs. (|10p . (jlip and (|12|) it can be inferred with mi = m 2 = m that the threshold 
behaviors of pp(s,mi,m 2 ) oc (s — 4m 2 ) 1 / 2 and of ps(s,mi,m 2 ) oc (s — 4m 2 ) 3 / 2 correspond to 
those of an S and of a P wave, respectively. Taking into account the intrinsic negative parity 
of the qq state, it implies the correct behavior under parity transformation of the bound state 
described by the qq state and its associated vertex (see Eq. (fTSj) ), 



III. STRUCTURE OF THE BOUND STATE FOR A SCALAR PARTICLE 

Assuming that the /o(980) scalar meson is made of components uu,dd and ss, one can 
decompose the total wave function as follows: 

^/o = -^={uu + dd) sin^mix + sscos6> mix , (13) 
V2 

or 

= * M sin # mix + * As) cos 9 mix = N s ( 4> {n) sin # mix + c/> (s) cos # mix ) , (14) 

where # m i x is the mixing angle between the nonstrange, ^ j n ) , and strange, ^ („) , flavor content 

Jo Jo 

of the wave function 2 . In what follows, unless otherwise stated, m n will denote the up or down 
quark mass (m n = m u = m<i), m s that of the strange quark and N$ is the normalization constant 
of the full wave function. 



A. The scalar particle on the light front 

The explicit covariance of this approach allows to write the general structure of the two-body 
bound state. For a scalar particle (see Fig. [T|) composed of a quark-antiquark pair of equal mass 



2 Consequently this implies a strange component for the wave function of the a, ty a = (uu + dd) cos# m i x /\/2 — 
sssin# m j x . However such a strangeness content does not seem to have an experimental support (see for instance 
Ref. This certainly points to a more involved structure of the a or /o(600) than that of a simple qq state. 



m q and four-momenta &2 and fei, we have (q = n or s) 

^ = i«ft)#(i,RiHy, (15) 
V2 

where w(A;i) and u(&2) are the usual antiparticle and particle Dirac spinors, and A^ q \x,Il 2 L ) is 
the scalar component of the wave function. Note that the color factor is not included in the wave 
function Eq. (|15p . Since the quark masses m q , in each component A^ q '(x, Rj_), are identical, the 
corresponding reduced mass is m q /2 and we chose the following Gaussian expression: 

A^\x,R 2 ± ) =exp(-16 v q k 2 q /mf), (16) 

where v q is a size parameter to be determined from experimental data and theoretical 
assumptions while the momentum squared, k 2 , given in Eq. © now reduces to 

Rj+m^x-l) 2 
k "" 4x(l-x) ' 1 j 

B. The scalar particle in the dispersion approach 

The soft constituent-quark structure of the scalar meson is given in this approach by the 
vertex 

«™^i> G „ (5) , (18) 

where Q a (— ^2) and Q a (k\) are the constituent spinor states of color a normalized by the color 
factor Nc = 3. For a scalar meson made of a quark-antiquark pair of equal mass m q , the wave 
function cj)^(s) of Eq. dHJ) can be parametrized as 

4,(^ a) = JL-^!L= Wq (k q ), (19) 

where k q is the modulus of the quark momentum in the center of mass momentum such that 

4k 2 q = s-4m 2 q . (20) 

The functional form (|19|) is so chosen as will be seen later, so as to simplify the normalization 
condition in Eq. Q. The function w(k) is defined to have the same functional expression as in 
CLFD: 

w q (k) = exp (-16 u q k 2 q /m 2 q ) , (21) 

and here again the size parameter u q is to be determined from experimental and theoretical 
considerations. 



9 



IV. PHYSICAL CONSTRAINTS FOR A NEUTRAL SCALAR 

As described in detail previously, the qq bound states are described in both formalisms by 
vertex functions which are related to Gaussian wave functions. These have to be normalized 
and their phenomenological size parameters determined. A standard approach chosen in the 
quark model is to calculate the decay constant with the appropriate loop diagram and fix the 
size parameter that enters the calculation so as to reproduce the experimental value of that 
constant (see Appendix |A|) . In this work, this is done for the pseudoscalar D- and 5-meson 
wave functions. However, the lack of knowledge of the experimental /o(980) decay constant 
makes it difficult to proceed similarly for the scalar meson. Furthermore, the mixing angle 
# m ix is not known a priori. We therefore resort to a different parametrization prescription by 
making use of D decay branching ratios which contain the /o(980) in the final state. In this 
section, we discuss the constraints on the scalar wave functions given by the normalization and 
the experimental data set chosen to determine the mixing angle # m ; x as well as the various size 
parameters in both formalisms. 

A. Normalization in CLFD 

According to the spirit of the constituent quark model, the state vector is decomposed into 
Fock components, and only the two-body component is retained. Since the state vector is 
normalized as 

(p',\'\p,\) = 2poS {3) (P-P')5 X '\ (22) 
it gives for a zero total angular momentum state the following normalization condition 45(: 



1=/. D(x,9,R±) ^ *A?A a *A*Aa' ( 23 ) 



'(*.«.R±) AiAa 
where D(x,9, Rj_), is the invariant phase space element given by: 



10 



Using the condition of normalization for the Dirac spinors, Y2\ u a(^) u b(^) = + m )ab an d 
J2\ v a(k)Vb{k) = — 7T '-)afe! we sum over all spin and color states and get for a qq component: 

E E <a 2 <i = ~ E E ^{kM^v^^ikM^ik^ 



color Ai 



color Ai,A2 



a,b,c,d Ai,A2 

\ (h + rn) da (A^) ab (h-m) bc (A (c,) U 



(25) 



= ~Tr [(fe + m)(^-m)(A(«)) 2 ] , 

where j4^^ is given by Eq. f)16[) . The result is similar for both the nn and ss components. 
There is no mixing term between the two components. With the scalar wave function written 
in Eq. (1141) . the normalization condition is therefore: 



1 = N$ 



ki [AM (x, Ri)) 2 sin 2 9 mix + (A^(x, R 2 )) 2 cos 2 9 



4mH 



4m? 



£>(x,0, R_l), 
(26) 



with k q given by Eq. (fT7|) . 



B. Normalization in the Dispersion Relation approach 



In the DR approach, for the scalar meson, the wave function [see in Eqs. ©, (|14p and (|19p ] 
is normalized according to Eq. ([S]). Taking into account the quark-content assumption of the 
/o(980) introduced in Eq. (|14p and making use of the form for cj)^ (s) or (f)^ (s) given by Eq. (|19p . 
the normalization condition for \I/j reads 



2 " +w 2 s (k) cos 2 8 mix ]dk, (27) 



l = iVf / fc 2 [u£(fc)sm „ mlx 

JO 

since the cross contributions vanish because of the orthogonality of the flavor states. In Eq. ([2j 
w n (k) is the nonstrange Gaussian component of <^ n '(s) whereas w s {k) is the strange one of 
4>( s \s), which implies two different size parameters v n and v s . The form (|2ip for ui q {k) leads to 
the normalization 

2 



iV 5 



7T 



1/4 



2 N§ 



32i/„ 



sm 



+ 



m 



9 \ — 
2 \ 2 



32i/ s 



cos 



(28) 



C. £) meson branching ratios 



The wave function for scalar meson /q(980) - denoted hereafter for simplicity /o - is 
constrained by the experimental branching ratios for the channels D + — > foir + ,D° — > 



11 



f K°,D+ -> f K+,D+ - f ir+,D+ 
are provided by different collaborations, E791 



BABAR [20], FOCUS [lj and E687 [2J,l2j: 



fpK + and .D — ► /ovr°. The experimental ratios 
23], ARGUS Q, 



21 



22, 



CLEO 



13, |4fl 



10, 



£791 [22, 23] : 



B(D + -» / vr + ) x B(/ -» tt+O =(1.9 ± 0.5) x 10" 4 , 



(29) 



ARGUS [19] : 
CLEO [17] : 
BABAR [20] : 



£(L>° -» f K°) x B(/ -> tt+tT) =(3.2 ± 0.9) x HT 3 , 



B(D - / ^ u ) x B(/ - vr+vr-) =(2.5^) x 10" 3 , 
B(D° -► /o^°) x B(/o -> K+iO =(1.2 ± 0.9) x HT 3 , 



(30) 



FOCUS [18] : 
FOCUS [18] : 



£(£>+ -► /o^ + ) x B(/ -► K + K~) =(3.84 ± 0.92) x 10~ 5 , 
B{D + / K+) x B(/ -► tt+vt") =(6.12 ± 3.65) x 10~ 5 , (31) 



£687 [24, 25] : 
£791 [21] : 
FOCUS [18] : 
FOCUS [18,29] 



B(D+ - /o7T+) x B(/ - 
B(Z>+ -> / vr+) x B(/ 
B(£>+ -> /Q7T+) x B(f 

B(D+ -> /o7T+) x fi(/o - 



K + K~) =(4.9 ±2.3) x 10~ 3 , 
-> vr+vr") =(5.7 ± 1.7) x 10~ 3 , 
-> Tr+vr") =(9.5 ±2.7) x 10~ 3 , 

K + K~) =(7.0 ± 1.9) x 10" 3 , 



(32) 



FOCUS [18] : 



£(£>+ - / K + ) x B(/ - K+K-) =(2.8 ± 1.3) x 10" 



(33) 



CLEO [46,47] : 



B(D° -► / vr°) ~ Bp vr+TT-vr ) x ^(D -► / tt ) = 
(1.1 ± 0.4) x 10~ 2 x (1.0 ± 0.8) x 10" 4 , (34) 



In Eq. ((31) F(D° -> / 7r°) represents the fit fraction of the (D° -> / vr ) decay [46j 3 . The 
/o(980) width is dominated by the /o decay into tttt and Combining their partial wave 



analysis 



48J of 



Collaboration 



2xo 



7r + 7r K + K with their study [491 ] of Xc ~> /o/o — * vr + 7r -7r + 7r the BES 



481 ] has determined the following ratio between the partial widths of the /o 



R 



r(/ -> vrvr) 



r(/ - ^vr) + r(/ - ifif) 



0.75 



-0.11 
-0.13- 



(35) 



3 The value used for B(D° — > 7r + 7r 7r°) is taken from Ref. [47] in order to be consistent with that of Ref. [46 
The more recent and precise value from Ref. [i| does not modify our conclusions. 



12 



Applying isospin relations, one finds the following branching fractions 

B(/ - tt+O = ^ = 0.50±°;^, (36) 

B(f - K+K-) = ^ = 0.1251H1. (37) 

The two-body branching ratios B(D — > foP) entering Eqs. (|29p to (|34p are then deduced from 
the branching fractions B(fo — > 7r + 7r~) and £>(/b — ► i^ + -ftT~) given in Eqs. (f36|) and (J371) . It is 
worth emphasizing that the results are strongly dependent on these branching fractions. Note 
that their experimental uncertainties are large. 



V. ELECTROWEAK AMPLITUDE 



In any phenomenological treatment of the weak decays of hadrons, the starting point is the weak 
effective Hamiltonian, which is obtained by integrating out the heavy fields from the standard 
model Lagrangian and reads 

H% C f =l = ^Y^VcKMCiMOifa) + h.c, (38) 

* i 

where Gf is the Fermi constant, Vckm contains products of the CKM matrix element, Ci(fi) are 
the Wilson coefficients, Oi(fi) are the operators entering the operator product expansion and \x 
represents the renormalization scale. In the present case, since we only take into account tree 
operators, the matrix elements of the Hamiltonian (|38p read, 

2 



{M l M 2 \H^f 1 \D) = ^Y, 



+ h.c, (39) 



V cq V: q Y J CimM 1 M 2 \Ol\D)^) 

q L i=l 

where q = d or s according to the transition c — > d or c — > s. The scale [i is chosen to be 
of order m c for D decays. The amplitudes (MiM2\Of\D) ([/,) are hadronic matrix elements, 
Mi and M 2 denote a pseudoscalar and a scalar meson in the final state. In Eq. (I39p the 
notation (M% M 2 \Of\ D)(n) reflects the fact that the hadronic matrix elements also depend on the 
renormalization scale fi. They describe the transition amplitude between initial and final states 
at scales lower than \x and give rise to the main uncertainties in the calculation, as they involve 
the nonperturbative regime of QCD. The operator product expansion divides the calculation of 
the amplitude A(D — ► M\M 2 ) oc Ci(^){MiM 2 \Oi\D)(p,) into two distinct physical regimes. One 
regime deals with hard or short-distance physics, represented by the Wilson coefficients Ci(fi) 
and calculated perturbatively, the other concerns soft or long-distance physics. The operators 
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Oi(n) can be understood as local operators which govern a given decay, reproducing the weak 



351 ] contain 



interaction of quarks in a point-like approximation. The Wilson coefficients Cj(/i) 
the physical contributions from scales higher than /i. Since QCD has the property of asymptotic 
freedom, they can be calculated in perturbation theory and include contributions from all heavy 
particles with m > fi, such as the top and beauty quarks, the bosons, and the charged Higgs 
boson. The dependence of the hadronic matrix elements and of the Cj(/i) on \i must cancel in 
the final decay amplitude which is a physical observable and thus scale independent. 

Working at tree level within the factorization formalism one obtains the following decay 
amplitudes: 

^VcdKd Kh - ^/„)«l/^ ° K+)sfn0 mix for D+^f 7r+, 
^V cd V: s (m 2 D+ -mj o )aif K ^ (m 2 K+ ) S m9 mix for D+^f K+, 



A{D -> f P) = { 



G 



2 VcsVJd ( m bo -mj o )a 2 f K ^ 



no . A n ) 



(m 2 KO )sm6 mix for D° -> f K°, 



%VcsV: d {ml t -ml)a 1 urf t ^ } »\ml + )co S e mh , for D+ -> 



^=V CS V* S (m 2 Dt - m 2 fo ) ai f K ^ s (m 2 K+ ) cos0 mix for D+ -> / ^ H 
li^cdKd Ko-m/ )«2W K o )sin0 mix for D° - / 7r°, 



(40) 



where aj(m c ) is written as dj for simplicity. In Eq. ()40p . and are the pion and kaon decay 
constants and 



oi(m c ) = Ci(m c ) + 



C 2 (m c ) 



a2(m c ) = C 2 (m c ) + 



Ci(m c 



(41) 



where iVc = 3. The flavor content u or s of the D and /q has been written explicitly in the scalar 
transition form factors (jrtp). With these factorized decay amplitudes, we can compute 

the decay rates using the following expression , 



T{D -> f P) 



8tt m 2 D 



\A(D -> / P) 



(42) 



where |p| is the modulus of the cm. momentum of the decay particles defined as 

y/[m* D - {mp + m f0 ) 2 ][m 2 D - (m P - m fo ) 2 ] 



IPI 



2m d 



(43) 



Finally, one defines the branching ratio B as the ratio between the decay rate T{D — * JqP) and 
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the total decay width Yd: 



B (44) 



A. Numerical inputs 

1. Values of CKM matrix elements and Wilson coefficients 
In the present numerical calculations, the CKM matrix elements, expressed in terms of the 
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5l| rely on the latest values extracted from charmless 



Wolfenstein parameters A and A 
semileptonic D decays Q|: 

A = 0.2257, A = 0.814. (45) 

The Wilson coefficients, at the mass scale fi = m c , are Ci(m c ) = 1.3777 and C2(m c ) = —0.6941 
(see Ref. {35] ]) from which we infer 

ai(m c ) = 1.1463, a 2 (m c ) = -0.2349. (46) 
2. Quark masses 

We use the subsequent standard constituent quark masses to calculate the transition form 
factors within the quark model approximation. 

m u = m d = 0.350 GeV, m c = 1.620 GeV, 

m s = 0.510 GeV, m b = 4.920 GeV. (47) 



[3] 



For meson masses, the following values [3] are used: 

m B ± = 5.279 GeV, m Ba = 5.369 GeV, 

m D ± = 1.869 GeV, m Da = 1.968 GeV, 

m D o = 1.864 GeV, m K ± = 0.493 GeV, 

m K o = 0.497 GeV, m fo = 0.980 GeV, 

m w ± = 0.139 GeV, m w o = 0.135 GeV. (48) 

The pseudoscalar decay constants fp are defined as usual by 

<P(pi)|gi 7 'V<Z2|0) = -ifppl (49) 
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Pi being the momentum of the pseudoscalar meson and the numerical values we used are, 

f n = 132 MeV, f K = 160 MeV, f D = 222 MeV, 
f Ds = 274 MeV, f B = 180 MeV, f Bs = 259 MeV. (50) 

The Fermi constant is Gf = 1.166 391 x 10~ 5 GeV -2 [jj and the world average D life-time values: 



r D o = 0.410 ± 0.001 ps, r D ± = 1.040 ± 0.007 ps, r D ± = 0.490 ± 0.001 ps. 



yield the total D decay widths Yd = 1/td. 



(51) 



VI. WEAK DECAY FORM FACTORS FOR P -> S TRANSITIONS 
A. Standard form factor notation 

The decays of b and c quarks are given by the weak current (even though only the 
^7^7 ?6(c) term is relevant in our case), 



Jt 



6(c 



) =9^(1 -7 5 )9 6(c); 



(52) 



where q is a light u, d or s quark. As usual, one can define the physical amplitude for a semi- 
leptonic decay X —*YIt/i by the expression 



M 



^i(S\J»\P)Jp, 



(53) 



where J^ cp is the leptonic current. In Eq. (|53p . (S\J^\P) is the hadronic matrix element including 
the weak current as defined previously. Introducing the total four-momentum K = P\ + P% and 
the four-momentum transfer q = P\ — Pi where Pi is the four-momentum of the pseudoscalar 
meson and P% that of the scalar meson in the final state, the hadronic matrix element can be 
decomposed as: 



{S(P 2 )\J»\P{P 1 ))=K»f + { q 2 )+q»f_(q 2 



(54) 



where f+(q 2 ) and f-(q 2 ) are the transition form factors and Pi and P2 are respectively the four- 
momentum related to the initial and final particle states of the hadronic current. Introducing 
then the scalar Toil 2 ) an d vector ^{q 2 ) form factors, the amplitude can be expressed as 



(55) 



(S(P 2 )|J^|P(Pi))=^i(g 2 ) 




+ Mi 2 ) 


\ K 2 " A 




q* 




q* 
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since K • q = M 2 — Mf , and M\ and M 2 being the masses of the initial and final meson. It is 
straightforward to derive the relationship between the two sets of form factors. One obtains 

Fi(q 2 ) = U(q 2 ), (56) 

F G tf) = f + {q 2 ) + -f-f-{q 2 ). (57) 
A • q 

Note that at q 2 = 0, JF^O) = T (0) = /+(0). 



B. CLFD formalism 



In the covariant light-front dynamics formalism, the exact transition amplitude does not 
depend on the light front orientation but in any approximate computation the dependence is 
explicit. However one can parametrize this dependence since the formalism is covariant. Hence, 
the approximate amplitude expressed in CLFD is given by the following hadronic matrix, 

(S(P 2 )\J^\P(P l )) CLFD = K»f + (q 2 )+q»Uq 2 )+^B(q 2 ), (58) 

where B(q 2 ) is a nonphysical form factor which has to be zero in any exact calculation. The last 
term represents the explicit dependence of the amplitude on the light front orientation m with 
ijj 2 = 0. In order to extract the physical form factor f±(q 2 ), without any dependence on u, from 
the amplitude (S(P 2 )\J^\P(Pi)) CLFD , we will proceed as follow. First, we calculate the scalar 
products X, y and Z which are defined by, 

X = K„. ( < S(P 2 )|^|P(P 1 )) CLirD = K 2 f + (q 2 ) + K ■ qf_(q 2 ) + K ■ uB(q% (59) 

y = q, ■ (S(P 2 )\J»\P(P 1 )) CLFD = K ■ qf + (q 2 ) + q 2 f„{q 2 ) + q ■ uB{q 2 ), (60) 

and finally, 

u-P l Z = ^-(P 2 | JHPi) CiFD = K ■ uf + {q 2 ) + q ■ uf^(q 2 ). (61) 
We define a variable y as the ratio between the scalar product of go ■ P 2 and u-Pi, 

• U "- /W - with P r P 2 = i(M 1 2 + M 2 - q 2 ). (62) 



y u-Pi M 2 + P l -P 2 ' 2 

Since P\ = (K + q)/2 and P 2 = (K — q)/2, we may also write 

_ ujJ^ _ u-(K-q) _ 4M| + K 2 - q 2 
V ~ ITFx ~ uj-(K + q) ~ AM 2 +K 2 -q 2 ' 



(63) 
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with uj ■ K = (1 + y)u:Pi and u ■ q = (1 — y)ujP\ . For q 2 > 0, it is convenient to restrict ourselves 
to the plane defined by u ■ q = 0. This condition is allowed in the system of reference where 
Pi + P2 = with P\q — P20 7^ 0- From the scalar products X , y and Z we can isolate the form 
factors f±(q 2 ) from B{q 2 ). Then, one gets the expressions for the form factors f±(q 2 )'- 



f±(q 2 ) = n*±(y,q 2 ,X,y,Z) 



where f2 is identical for both form factors f±(q 2 ) and can be written as, 







1 2 > 



{M 2 (y -l) + q 2)y-(y- l)M 2 \ [(y ~ l)K + (y + l)q 
where the functions ^f±(y, q 2 , X, y, Z) are: 

(y, q 2 ,X, y, Z) ={y 2 - l)X + (y + l) 2 y + [(1 - 3y)M 2 - (y - 3)M 2 2 + (y - l)q 2 
V + (y, q 2 ,X, y, Z) ={y - 1) 2 X + (y 2 - l)y + [(y - l)M 2 - (y - 1)M 2 2 + (y + l)q 2 

or in terms of the variables K and q, 



(64) 



(65) 



Z, 
Z, 



(66) 



*_(!/, (? 2 , X, y, Z) =(y 2 - l)X + (y + lfy+ [(1 - y)K 2 - (1 + y)K-q 
V + (y, q 2 ,X, y, Z) ={y - l) 2 X + (y 2 - l)y + \(y - l)K-q + q 2 (y + 1) 



Z, 
Z. 



(67) 



The second step is to express the amplitude {S{P2)\J^\P{Pi)) CLFD without using the form 



factors f±(q 2 )- The leading contribution to the transition amplitude (S(P2)\J^\P(Pi)} 



CLFD 



is 



given by the diagram shown in Fig. By using the CLFD rules, one can derive the matrix 
elements from the diagram (Fig. [2]) and one has, 



(S(P 2 )\J^\P(P 1 ))^ LFD 



D(x,6,-R ± )Tr 



V2 



As(x ',Rl 2 )(mi + jti)ij^, 5 (m 2 + fc) 



x ±A^'\x,R 2 ± )(m 3 -h) 



1-x' 



where A^ q ^(xjR^) and As(x' ,Tt 2 ) are the pseudoscalar and scalar wave functions defined in 
Eq. (|A2|) and Eq. (|16p respectively. Note that x and x' are the fraction of the momentum carried 
by a quark q% (the spectator quark) as given by: 



x = — — . and x = — — , then y 



x 



U)-P\ UJ-P 2 X' 

and one also has = Rj_ — x'q. Now, one can replace the hadronic matrix element 
{S(P 2 )\J^\P(Pi)) CLFD , which appears in the scalar products X,Y,Z (Eqs. §59\ [ffil EI])), 
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by the hadronic matrix elements (S(P 2 )\J fl \P(Pi)) < ^ LFD calculated by applying the CLFD 
diagrammatic rules and given in Eq. (|68p . Hence, by using Eq. (|64p we are able to compute the 
form factors f±(q 2 ) as a function of q 2 and this over the whole available four momentum range 

< q 2 < 9max- 



C. Dispersion relation approach 

The pseudoscalar to scalar transition amplitude is calculated from the triangular Feynman 
diagram shown in Fig. O where also the kinematical variables are displayed. For the evaluation 
of the space-like transition form factor (q 2 < 0) the internal constituent quarks are put on-mass 
shell. Moreover the external momenta are considered off-shell with 

Pf = Sl , Pi = s 2 , (Pi - P 2 f = q 2 . (70) 



To derive the transition amplitude (|54|) we need the constituent quark matrix element of the 
weak axial current which we write 

(QUh)\qi(0)(-h^ 5 )q2(0)\Q a 2 (k 2 )} = -ihx{q 2 )Ql{k x )^Q 2 {k 2 ). (71) 

The function f 2 i(q 2 ) is the constituent quark transition form factor. Since no formal derivation 
of the quark model from QCD exists, it is unknown. In the following we make the assumption 
/21 ~ 1 and drop the factor altogether owing to the fact that constituent quarks behave very 



much like bare Dirac particles 



52] 



In the DR approach the transition form factors f±(q 2 ) of Eq. (|54p are expressed through the 
double spectral representations: 



2 I ds 2 G V2 (s 2 ) f dsi G Vl (s\) . . 2 , ,-.->, 

f±(q ) = / — f / 772V A ±(*i> s 2,9 ;rn 1 ,m 2 ,m 3 ). (72) 

7r(s 2 - M|) J 7T(S1 - M{) 



The functions A±(si, s 2 , q 2 ; mi, m 2 , 1713) in the above equation are the double spectral densities 
of the triangle Feynman of Fig. [3] in the P 2 - and P 2 - channels. They can be obtained [42] from 
the following equation 

(Pi + P 2 ) fl A + (s 1 ,s 2 ,q 2 ;m 1 ,m 2 ,m 3 ) + (Pi - P 2 )^A^(s 1 ,s 2 ,q 2 ;m 1 ,m 2 ,m 3 ) 

= — f d^kid^k^hSikf - ml)5(k 2 - m 2 2 )5(k 2 - m|)5(Pi - k 2 - k 3 )5(P 2 - k 3 - h) 
8vr J 

x Tr [-(#! + m 1 ) 7 ^7 5 (^ 2 + m 2 )i 7 5 (m 3 - Jk 3 )i] , (73) 

where m 2 > mi. Explicit expressions for A±(si, s 2 , q 2 ;mi,m 2 ,m 3 ) are given in AppendixO An 
analytical continuation in q 2 allows us to write the transition form factors for q 2 < (m 2 — mi) 2 
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as 



oo s+(s 2 , q 2 ) 

2 f ds 2 G V2 (s 2 ) f dsi G Vl (si) B±(si,s 2 ,q 

f±{<1 ] = J *(s 2 -M 2 ) J 



(mi+m 3 ) 2 Sj (s 2 ,g 2 ) 



16n( Sl -Mf) \ 3 / 2 (si,s 2 ,q 2 ) 



9fj( 2, f ds 2 G V2 (S 2 ) f 

[Q ' J <s 2 -M 2 ) J Wn(si - sf )3/2 

4(q 2 ) s?(s 2 ,q 2 ) 

G Vl ( Sl ) B±{s u s 2 ,q 2 ) G Vl {sf) B±(s?,s 2 ,q 2 ) 



(74) 



.(Sl - s[) 3 /2( Sl _ Ml) (sf - sf) 3 / 2 (sf - M 2 )_ 

The functions sf(.S2, <? 2 ) = {■\fs 2 — \fq 2 ) 2 and Si(s 2 ,q 2 ) = (^fs^ + \/q 2 ) 2 are the roots of 
\{s\,s 2 ,q 2 ) = {s\ + s 2 — q 2 ) 2 — ^s 2 s\. The expressions for B±{s\, s 2 , q 2 ) are given in Appendix |B1 
along with the integration limits s^(s 2 , q 2 ), s~l(s 2 ,q 2 ) and s^iq 2 )- 

We note that although the diagrams for D — * /o(980) and ir — > /o(980) are very similar in 
the calculation of their spin trace, the main difference is that no kinematical factor is involved 
at the triangle apex where the interaction vertex 7^(1 — 75) in the case of the heavy meson 
decays stems from the weak interaction. In the tt — ► /o(980) form factor all constituent masses 
are identical, while in the present case there are two mass scales, namely the charm and a 
light or strange quark. It can be seen, in Appendix [HI that the expressions for the functions 



B+(s-\ , s?,, q 2 ) which enter the spectral densities vanish identical 



in Ref. 



531 ] . which relies on the method developped in Refs. 



54j, 



y for mi = m 2 = 7713. Indeed, 



551 ]. the transition amplitude 



for 7r — > /o(980) calculated on the light front was shown to vanish as t ~ q 2 ^ — * 0. This is in 
contradiction with experimental findings in ir~p — > 7r°7r°n reactions. We can make the parallel 
for the behavior of our transition amplitude in the limit q\ — * and confirm the vanishing of 
our form factor for t = q 2 — > if all internal quark masses are equal. Thus, had we calculated 
the transition to /o(980) from a pion, we would obtain f±(q 2 ) = for q 2 — > 0. We ascribe this 
discrepancy to our simplified qq picture of the /o(980) whereas other contributions, likely from 
pion and kaon clouds, may modify the form factors in particular at low momentum transfer. 
For m 2 7^ mi = 7713, however, we deduce from the expressions of our dispersive representation 
that the integrands in Eq. (|74p do not vanish for q 2 — ► 0, nor do the integrals as confirmed by 
our numerical calculations. 

As mentioned before, the form factor in the region < q 2 < (rn 2 — mi) 2 can be obtained 
by analytic continuation of the expression in Eq. (\72h for q 2 < 0. In this space-like region, 
the function A(si, s 2 , q 2 ; mi,m 2 ,ms) in Eq. (|B6|) has no square-root cuts related to the zeros of 
A 1//2 (si, s 2 ,q 2 ) (they lie on the unphysical sheet) and both form factors are given by just the first 
term in Eq. (f74l) . Note that the vertex functions G v (s) are not singular for s > (m 2 + 7TI3) 2 and 
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s > (rrii+mz) 2 and that B±(si, s%, q 2 ) are polynomials. Thus, the analytic properties of the form 
factors are determined by the sole behavior of the function A 1 / 2 (si, S2, q 2 ) for positive q 2 . One 
may study the structure of the singularities of the integrand in the complex s± plane for a fixed 
real value of S2 > (mi + m.3) 2 , which implies external S2 integration and internal s\ integration 
(interchanging the integration order leads to an equivalent integration contour). At q 2 > 0, the 
square-root cut endpoint sf 1 moves onto the physical sheet through the interval from to sf 
to the left of . This occurs for a value of S2 > -^(q 2 ), where s^iq 2 ) is obtained as the solution 
to the equation sf(s2,<? 2 ) = s±(s2,q 2 ) and given in Appendix [Bl The integration contour of si 
in the complex plane must be deformed so it encompasses the points s;p and . It therefore 
contains two integration segments, one being the normal part from to s~f and the other the 
anomalous part from to s7. The double spectral density for this anomalous part, on the other 
hand, is obtained from the discontinuity of the function A 1//2 (si, S2, q 2 ) which can be written as 



A 1 / 2 (si, S2, q 2 ) = J{s x - sf)(si — sf), bearing in mind that the branch point lies on the 
unphysical sheet. Hence, one has to calculate the discontinuity of l/(sx — sf ) 1//2 which is just 
twice the the function itself 5(| . This explains the integration limits and the factor two in front 
of the second integral in Eq. (|74p . The subtraction term in the third line of Eq. (|74p stems from 
the function 1/A(si, S2, q 2 ) that enters the complete expression for A±(si, S2, q 2 ] mi, rri2, m.3) in 
Eq. (IBlj) and which is singular in the lower integration limit sf~. It was shown in Ref. [57j] that 
an accurate application of the Cauchy theorem yields this subtraction term. 



VII. NUMERICAL RESULTS 
A. The fit procedure 

As we have discussed in the preceding sections our final aim is to predict form factors for 
Bfg\ — > fo transitions. To achieve this goal, we first have to acquire a good knowledge of the /o 
wave function. This will be done through the evaluation of theoretical branching ratios [Eq. (|44p ] 
for Dr s \ — > fo transitions, which implies the calculation of form factors that rely on the /o wave 
function, as can be seen in Eqs. (|68p for CLFD and (174ft for DR. Since on the one hand meson 
masses and decay constants are measured, and on the other hand constituent quark masses as 
well as Wilson coefficients are known from theoretical considerations and given in Sec. |Vj the 
evaluation of the branching ratios depends only on the /o wave function parameters: two size 
parameters u n , u s and the mixing angle m [ x . The overall normalization N$ in Eq. (114ft is fixed 
by means of Eq. (|26ft for CLFD and Eq. (|28ft for DR. Once the fo wave function parameters are 
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given, the form factors J- Q (s) °(q 2 ) and hence the branching ratios can be determined. These 
parameters will thus be constrained, via a least-square x 2 fit 4 ! by the experimental branching 
ratios given in Eqs. (|29p to (|34p . Note that there are two equivalent solutions for # m i x , as the 
mixing angle enters quadratically into the decay rate formula Eq. (|42|) . As an additional physical 
constraint, we choose to impose the relation 



between the strange and nonstrange components of the /o(980) wave function. This forces the 
strange component to be wider in momentum space, the size parameter v a being divided by m 2 
in the Gaussian wave functions given in Eq. (|16p or Eq. (|21|) . assuming that \ss) is more tightly 
bound and compact in configuration space. This effectively reduced parametrization proves to 
be decisively more stable, while not spoiling the fit. 

We will see that this first simple approach, a two parameter fit attempting to reproduce all 
data listed in Eqs. ([29]) to (f34"|) . to which we will refer to as fit 1, provides, partly because of 
the large experimental errors, a fair agreement with the data though not entirely satisfactory. 
Indeed, so far we a priori miss relevant physics in these decays such as corrections to simple 
tree-order topologies. We must include higher-order and power suppressed contributions in the 
appropriate channels. We here consider penguin and annihilation topologies which we now 
discuss in turn. 

In the decays discussed in Sec. IIVC1 penguin topologies only contribute to the D + — > /o7r + , 
Df — > foK + and D° — > fair amplitudes. The magnitude of the CKM matrix elements 
implies that for charmed penguins the penguin contributions can be and are usually discarded 
since VcdV* d ~ V CS V* S is three orders of magnitude larger than V c bV* b . Nonetheless, in order to 
try to even more constrain the scalar mixing angle, we have inserted phenomenological penguin 
amplitudes where they are operative. 

We have parametrized these contributions by a universal amplitude so that we have only 
modified the linear combination of Wilson coefficients, af. 



v s = v, 



(75) 




(76) 



4 The routine MINUIT [58| has been used to minimize the y 2 in this work. 
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which leads for the amplitude A(D + — > /ovr + ) to the substitution 



A(D + - / o7 r+) => .4(£>+ - /ovr- 



+ 



2 



KdKTd - m %)f^o U (""4+) sin 9mi x X p (p p , 6 P ), (77) 



and similarly for the other two channels with the same X p (p p ,5 p ) = p p exp(i(5 ; 



As has been argued in Ref. 
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vi- 



59( | . weak annihilation amplitudes are not negligible for 
the decays D — > PP,SP and are comparable to the tree amplitudes. This occurs because 
these annihilation amplitudes, denoted in the literature by W exchange or W annihilation 
topologies, can receive contributions from long-distance final-state interactions. At the hadronic 
level, the quark rescattering is manifest in s channel resonances and the W-exchange topologies 
receive contributions from, for example, the 0~ resonance X(1830) [59]. Thus, we introduce 
a phenomenological annihilation term, X a (p a ,S a ), in the D° -» f (980)K° decay channel such 
that 

A(D° -» / (980)K°) A(D° /o(980)A'°) + G F X a ( Pa ,S a ) (78) 

with X a (p a ,5 a ) = p a exp(i^ a ). The modulus p a and phase 6 a are free parameters, the natural 
scale of p a is in principle given by the decay constants fp, f% and // . We stress that neither the 
contribution from penguin nor from annihilation amplitudes will allow to resolve the ambiguity 
on the mixing angle m [ x . 

In the following we introduce the effective transition form factors which, in the nonstrange 
sector read, 

<r /o (« 2 )=<r /o (^)^ (79) 



and in the strange one 



<f J V) = ^o P f /( V) cos0 mix . (80) 



In principle, the six parameters should be fit to the branching ratios listed in Eqs. (|29l) to 
(I34h . It turns out, as expected from the arguments given above, that in both approaches, CLFD 
and DR, the contributions of the penguin amplitudes are vanishingly small and do not lead to 
any improvement of the fit while the mixing angle maximally changes by 1°. In fact, the phase 
of the penguin amplitude is nearly zero and the modulus is very small. We conclude that we 
may just ignore its contribution. We will therefore refer from now on to fit 2 as a four parameter 
fit which includes solely the annihilation amplitudes as correction to the tree level. 
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Before discussing in details the results of our calculations, we wish to point out the large 
experimental errors that appear in the constraining data. There are furthermore inconsistencies 
in these data as can be seen for instance in the FOCUS experiment [18|], for the — > /o7r + 
channel. Here, we observe a discrepancy in the decay magnitude between the channels where 
the Jo (980) decays into a two-pion or two-kaon pair as well as in their errors. Partly, this may 
be ascribed to the use of the different branching fractions B(/o - > 7r + 7r~) and B(fo — * K + K~). 

Considering the theoretical ratio Ri = A(D + — > foTr + )/A(D + — > /qK + ) and the 
corresponding one for the Df meson, R 2 = A(Df — > foir + )/A(Df — ► JqK + ), one observes 
that they are equivalent when working at the tree level approximation for the decay amplitude 
if one assumes that F^ h {q 2 ) has roughly the same value for q 2 = m 2 and m 2 K and similarly 
for J~n a ~* {(j 2 )- Using Eq. (|40p for the decay amplitudes of these channels, the ratios R\ and 
i?2 are proportional to the same CKM matrix elements, V* d and V^ s , and to the pion and kaon 
decay constants; they are of the order of 4. 

Experimentally, though, this order of magnitude is strongly violated when data from FOCUS 
(BTl{D + -» f K+) = (3.07±1.65)xl(r 4 ), fromE687 {B1l(D+ -> / vr+) = (3.92±2.63) x 10~ 2 ), 
as well as from FOCUS (BK(D+ / 7r+) = (5.60 ± 3.08) x 10" 2 ) are used. These data appear 
to be incompatible with the other data. Hence, we shall study two cases in the four parameter 
minimization space, one with 12 data, referred to as fit 2a, the other with 9 consistent data 
referred to as fit 2b. 



B. The / wave function 

Table U which corresponds to fit 1 (with 2 parameters and 12 branching ratios) shows that 
the factorization model at the tree level order allows for a fair representation of the data with 
reasonably well defined parameters v n and # m i x given in Table ITU An obvious discrepancy occurs 
for the Dq — ► foK° channel, the apparent agreement with the BABAR data being only due 
to the very large experimental error. 

The stability of our fit is illustrated in Fig. [4] for both approaches. The x 2 /d.o.f. function 
is, in both cases (CLFD and DR), smooth and has well defined minima as a function of the 
mixing angle 9 m i x . We find a mixing angle # m ; x = 32° ± 4.8° with the CLFD model and 
#mix = 41.3° ±5.5° for the DR model and the symmetric angles with respect to 90°. These va' 
are in fair agreement with the ones estimated from Df — > /o(980)-7r + and Df — > (frir + decays 3|] 
which cover the rather wide range 20° < 9 < 40° and 140° < 9 < 160°. 

The addition of an annihilation amplitude in that channel D° — ► /o(980)-RT° does considerably 



ue 

n 
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improve the quality of the agreement with the complete set of data (fit 2a) , as seen in Table IIIIl 
The results for the parameters u n and # m i x (Table HV| are extremely stable as compared to those 
of fit 1 (Table P). 

Finally, retaining only the nine consistent data as explained above, we obtain (fit 2b) a further 
improvement of the x 2 /d.o.f. illustrated in Table PVT. The wave- function parameters (Table IVT1) 
for the CLFD model are stable as compared to those in Tables [XT] and IIVI whereas for the DR 
model, the range parameter increases by about 20% while the mixing angle remains stable. As 
for fit 1, the stability of fit 2b is illustrated in Fig. 0m both approaches. The x 2 /d.o.f. function 
is, in both the CLFD and DR models, smooth and has well defined minima as a function of the 
mixing angle 9 m \ x . 

The prediction for # m i x differs by about 10° in the two approaches. This can be explained as 
follows; in both approaches we employ equal Gaussian parametrizations of the vertex functions 
introduced in Eqs. (|16p and (|2ip . yet the dynamics that enters the loop diagram associated with 
the meson normalization differs somewhat in each case. In particular, in the DR approach the 
condition in Eq. (|9|) implies a vertex renormalization due to soft rescattering of the constituent 
quark in the vicinity of the meson pole mass M 2 . These modifications in the calculation of the 
normalization already cause differing normalization values in the case of the heavy pseudoscalars. 
As seen in Table ITXl in Appendix El although the values of v are very close in both models, the 
normalizations are quite different. This feature of the normalization is even more apparent for 
the /o(980) but, in addition, the size parameters v are an order of magnitude apart which results 
in different weights and ranges of the bound state vertex functions. These unequal weights enter 
the form factor calculations where they are compensated by the different normalizations N$ in 
both models. However, another degree of freedom comes into play here, namely the /o(980) 
mixing angle whose value can also compensate the Gaussian weights and thus competes with 
the normalization. Since for small momentum transfers, q 2 = m 2 and m 2 K , the effective form 
factors of CLFD and DR [see Eqs. (|79p and (|80p ] must be very close in order to fit the data, 
the product of Ns, the Gaussian weights and the sine or cosine of the mixing angle in the decay 
amplitudes Eq. (|40p must agree up to small variations inherent to a fit with two different models. 
The normalization Ns being not equal in CLFD and DR, this results in the observed variation 
of about 10° in the mixing angle. 

All fits of the branching ratios only constrain the /o(980) wave function at very small relative 
momenta k 2 , of the quark pair as given in Eq. (|20p . Though the introduction of the annihilation 
amplitude considerably improves the fit, its consequences on the scalar meson parameters are 
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rather limited. 

C. P — ► S transition form factors 

With the parametrization of the scalar-meson wave function in Table IVu resulting from fit 2b, 
we compute the pseudoscalar to scalar transition form factors D — * /o(980)^ n \ D s — * /o(980)W 
and can now predict B -» / (980)( n ) and B s -> / (980)^ . Indeed, with the values of Table [IX] in 
Appendix [Aj we can compute, employing Eqs. (|56|) . ([57]) . ([64"]) and ([73]) . the PS 1 — > S transition 
form factors D -> / (980)( n ), D s / (980)( s ), 5 -> / (980)( n ) and B s -> / (980)( s ) for any 
kinematically allowed momentum transfer q 2 . In the CLFD formalism this is done for q 2 > 
whereas in DR these form factors are evaluated for spacelike and timelike values of q 2 . Since 
we compare the two models, we only consider the positive range of q 2 . The momentum-transfer 
dependence of the effective Fo(q 2 ) and Fi(q 2 ) form factors is plotted in Figs. [6] and [7J for the 
D and D s transition form factors in both models. It is worthwhile to mention that, in the DR 
formalism, the anomalous contribution to the form factors, namely the second term in Eq. (|74|) . 
only sets in for momenta of q 2 > 0.6 GeV 2 in Dr s \ — > /o(980) transitions. Therefore, in the 
momentum range of interest here, this contribution is negligibly small and only the Landau part 
of the integrals is of interest. To make these effective form factors readily available, we assemble 
in Tables IVTI1 and IVIIII a list of their values for a few specific values of q 2 , namely q 2 = m 2 , m 2 K 
and m 2 . 

As can be read from Table IVIII and Fig. and [7J both models are in fair agreement for 
the range of timelike momenta < q 2 < 0.1 GeV 2 in the transitions D — ► /q(980)^ and 
Ds -> /o(980)( s ). This is expected as in the fit we fix the model parameters via the effective 
form factor Fo(q 2 ) for q 2 = m 2 and m 2 K barring any other changes in the decay amplitudes of 
Eq. (|40l). 

For the B to scalar transitions 5 , the kinematically allowed range is much larger than extending 
the momentum transfer squared up to 15 GeV 2 . Hence, once again, we do not consider 
contributions of the anomalous term in Eq. (|74p in the DR formalism. The effective form factors 
F (q 2 ) and F^q 2 ) are plotted in Figs. [8] and H Table EHD gives a few values at q 2 = m 2 , m 2 K , m 2 
and m 2 D . In Fig. [HJ for the B — > /o transition, one observes similar results to those obtained for 
the D to scalar transitions, whereas for the B s — * Jq, the difference between the DR and CLFD 

5 Our attention has been drawn by R. Dutta 60] to a work with S. Gardner where they obtained similar results 
with the use of the constituent quark model combining heavy quark effective theory with chiral symmetry in 
the light quark sector. 
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predictions is considerable as can be seen in Fig. [9l 

The magnitude of the slopes for Fo(q 2 ) and F\{q 2 ) point at different dynamical features for 
larger q 2 despite the use of similar vertex functions in both CLFD and DR. This is true in 
particular for large q 2 ~ m 2 values in B — > /o(980) transitions where one expects perturbative 
QCD effects to be relevant. It is likely that the Gaussian vertex form of the Bethe-Salpeter 
amplitudes which describe both the heavy pseudoscalar and the light(er) scalar bound states 
are not appropriate at large momentum transfers. In the D decays, the differences are even 
more pronounced — whereas at the maximum recoil point q 2 = the DR approach values for 
Fi(q 2 ) are slightly larger in magnitude than those from CLFD, they evolve more slowly and at 
(/max the CLFD predictions are considerably larger, as seen in Figs. [8] and [9j In this case, the 
momentum transfer range < q 2 < 0.6 GeV 2 is lower than the meson mass m 2 D and the process 
should be more dominated by soft physics. Therefore, the deviations between DR and CLFD 
cannot be ascribed to the behavior of the vertex functions and are intrinsic to the dynamical 
assumptions in either model. A feature of the DR model is that the function f-(q 2 ) decreases 
more rapidly than f+(q 2 ) increases, in particular for the D — > /o(980) transition form factors. 
This steeper slope as well as the factor q 2 /{m 2 D — m 2 o ) which is larger than q 2 /{m 2 B — rn 2 Q ) in 
Eq. (|57p also explain the negative slope of Fo(q 2 ) for the D — > /o(980) transitions. However, the 
difference with the CLFD form factor prediction is striking and only in the momentum domain 
of the pion and kaon mass can agreement be found. The problem of model dependence appears 
at larger momentum transfer, where various models yield rather different results, whereas at 
q 2 = Ref. 60(] seems to confirm our results. 

Regarding the general behavior of the transition form factors, in DR one observes that 
they are very sensitive to the function b± which strongly depend on the quark mass difference 
(Eqs. (|B4j) and (1B5I) ). In CLFD, the form factors are controlled by the function Q, introduced 
in Eq. (|65p . which forces Fo(q 2 ) and F\(q 2 ) to behave as l/(a + flq 2 ) and therefore become very 
large at the kinematical limit whenever the denominator tends to zero. 

It is worthwhile to recall that quark model predictions have a constituent mass dependence 
causing a systematic error in the computation of the form factors. This is in particularly true 
for the light sector where it is known that the dressed-quark mass receives strong momentum- 
dependent corrections at infrared momenta, an expression of dynamical chiral symmetry 
breaking. The enhancement of the mass function in the light-quark propagators is central to the 
occurrence of a constituent-quark mass scale. On the other hand, the impact on heavy-quark 
propagators of chiral symmetry breaking is much less marked for c-quarks and even less so for 
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6-quarks. It can be shown that that the heavy propagator S(p) = (]/> — ttiq) is justified for b 



quarks and to a certain extent also for c-quarks 



6 1| . Thus, in the approach of same propagators 



for light and heavy quarks, with a light constituent mass of m u <j = 0.35 GeV, a certain mass 
dependent uncertainty is implicit. We also remind that both D and B mesons are lightly bound 
and that the bound state condition M 2 < (m± + m-2) 2 is only fulfilled in the quark model if the 
light-quark mass is chosen to be large. However, since we make use of the features of confining 
models, this constraint does not affect our predictions. 

As an example, if we choose for the light-quark mass m U ( ^ = 0.25 GeV, modifications of 
the form factors magnitude at larger q 2 values are not insignificant. In the DR approach, 
for instance, a decrease of the light-quark mass, which implies a readjustment of the meson 
parameters to fit their decay constants, Fo(q 2 ) and F\{q 2 ) evolve more rapidly and overall we 
observe modifications of the order of 10% for q 2 up to the squared kaon mass. Changes in the 
strange quark mass scarcely alter these form factors on the other hand. This observation is more 
striking for D — > /o(980) transitions, where the heavy-light quark mass differences (m c — m u ) 2 
and (m c — m s ) 2 are smaller than when a 6-quark is involved. A proper treatment of dressed 
light-quark propagators should remedy this situation. 



VIII. EPILOGUE 



We have investigated the role of the scalar meson /o(980) in quasi- two-body decays of Di s \ 
and Bt s \ mesons focussing on the weak transition form factors Dr s \ — > /o(980)^^ and Bt s \ — > 
/o(980)^\ which are of particular interest to flavor physics. In order to obtain a consistent 
parametrization of the /o(980) wave function, we first applied a simple factorization ansatz to 
these D decays where the approach is reasonable. Here, the quasi-two body — > /o(980)P 
branching ratios are deduced from the experimental ones for — > tttttt, KKir and the 
knowledge of the /o(980) — > n + ir~ and /o(980) — > K + K~ branching fractions. Once the 
scalar-meson parameters are determined by fitting the matrix element — > /o(980)( s ^ to 
experimental data, they are readily available for other flavor changing matrix elements involving 
the 6-quark although in that case the approach is on less firm grounds. The short-distance 
physics in the factorization is known from perturbation theory applied to the operator product 
expansion and codified in terms of Wilson coefficients. The long-distance effects concern two 
sets of form factors; namely the experimentally known decay constants and the heavy-to-light 
transition form factors. The latter are nontrivial objects which involve quark as well as hadron 
degrees of freedom. In our approach, we have modeled these form factors with triangle diagrams 
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(at the tree level) in the impulse approximation. The mesonic Bethe-Salpeter amplitudes are 
described by Gaussian two-quark vertex functions which introduce size parameters. In the case 
of the scalar meson, we also need a mixing angle between the strange and nonstrange components 
of its wave function for which we assume the simplest possible quark structure. That is to say, 
we neglect higher Fock states or possible hadronic dressings which may enrich the qq state with 
other components such as \KK), \ivrf) etc. in order to perform an actual calculation. As noted 
previously in Sees. I and II, a consequence of the mixing is the presence of a strange component 
in the a or /o(600) state, strange content which does not seem to be experimentally observed. 
A specific discussion of the structure of this broad state is outside the scope of the present study 
and would require, as we just pointed out, to work beyond the simplest two-quark structure. 

In this work, we have examined two different but explicitly covariant approaches to 
establish the model dependance of the form factors. In both model calculations, the impulse 
approximation is used and quark masses as well as dynamical assumptions are the identical, 
though certain kinematical aspects differ. In particular, in the DR approach internal quarks are 
put on-mass shell and the amplitudes are expressed as double dispersive integrals of the triangle 
diagram's discontinuity over initial and final mass variables. In contrast to the DR approach, in 
the CLFD calculation the integration is performed over the internal loop momenta. Moreover, 
even though the Bethe-Salpeter amplitudes of the and /o(980) have identical Gaussian 
forms, the meson vertex normalization is not identical in both models. 

These differences may be the origin for certain discrepancies we find in our results. In 
fitting the set of experimental — > /q(980)P branching ratios, we do obtain similar values 
for the mixing angle. Overall, the fit quality is comparable and rather good given the large 
experimental errors. However, while at small momentum transfer, around the light meson 
masses m 2 and m 2 K , we find very similar transition form factors, for larger values of q 2 where 
no experimental constraints exist the discrepancy is obvious. In the case of Br s \ — ► /o(980) 
transitions, stronger deviations between both models are observed. For the Dr s \ — > /q(980) 
transitions, the discrepancy is already obvious for q 2 < m 2 K as seen, in particular, in the different 
slopes of Fo(q 2 ) obtained in DR and CLFD. This is also a hint that the constituent quark model 
may be reliable solely for a certain domain of q 2 . 

Clearly, dynamical aspects of QCD, such as running quark masses, are important in the 
computation of these form factors. 

Furthermore, the parametrization of the heavy mesons depends on the precise knowledge of 
the pseudoscalar decay constant. As confinement is only approximately achieved and dynamical 
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chiral symmetry breaking not realized in either model calculation, some of the uncertainty defies 
any quantification. When these formalisms are applied to calculations which can be compared to 
observables such as decay constants, typical deviations from the experimental values are of the 
order of 10% — 15%. Given the large errors in the experimental D( s ) — ► tttttt, KKtt branching 
fractions and the still elusive structure of the scalar /o(980), assuming a qq composition, this 
provides a lower bound of our theoretical error 6 which we estimate to be of the order of 25%. 

Nonetheless, we consider that there is a domain of validity for these models which overlaps 
with the typical momentum transfers q 2 that occur in leptonic as well as nonleptonic weak decays 
of -D( s ) and mesons. The present study provides a first calculation of heavy pseudoscalar to 
scalar meson transition form factors at the exact momentum transfer values q 2 = m 2 ^, m 2 K , m 2 
and mfj without resorting to any extrapolation. Surely, this work leaves plenty of room for 
improvement; obviously a better understanding of the scalar-meson structure is of foremost 
concern, but a more genuine realization of confinement and dynamical chiral symmetry breaking 
is also desirable. 
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6 which can be roughly tested by varying the size parameters and mixing angle within the error ranges shown in 
Table [VI] 
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APPENDIX A: PSEUDOSCALAR MESONS IN THE QUARK MODEL 



1. CLFD 



For a pseudoscalar particle composed of an antiquark and a quark, of mass mi and m 2 
respectively, the general structure of the two-body bound state has the form: 

1 



rj/f ) = N P ^ {qq ] 



. will! > = —u(k 2 )A^ ; (x,Ri) 75 v(h), 



|(99 )f„ T?2 



(Al) 



where v{k\) and u(k 2 ) are the usual Dirac spinors, and Ap (x, Rj_) is the scalar component of 



the wave function written as 

|(99 )(„ -r>2 \ _ A(m)(iJ2\ 



A™ >(x,Ki) = A^ q >(k') = exp(-4u^/m( 2 ), 



(A2) 



where Np and v are parameters to be determined by comparison with experimental data; the 
reduced mass is m\ 2 = m\ra 2 j{rn\ + 771,2) and k 2 is given by Eq. ([6]). For the pseudoscalar 
mesons we make use of the experimentally well established values for their decay constant. 

In CLFD the normalization condition for a pseudo scalar meson of zero total angular 
momentum reads as follows : 



1 



/ . D(x, 9, Rj_) J2 ^TaM^*' 
j(x,e,R ± ) A1A2 



where, in close analogy with Eqs. (|25p . one has, 



Am) = ^p Ti 

/ j ' Ai,A2 t / Ai,A2 2 
Ai,A2 

so that, finally, 



l = Np I D(x,e,Tlj 
l(x,§,-R ± ) 



(99 )/ 



R 2 + (xm 2 + (1 — x)mif 
x(l — x) 



[4 99) (x,r; 



2 \i 2 



(A3) 



(A4) 



(A5) 



where one recalls that D(x,6,H±) is the invariant phase space element already defined in 
Eq. tm. 



2. Dispersion approach 



Similarly, the the two-body bound state for pseudoscalar meson is given here by 



(A6) 
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where Q a (ki,mi) represents the spinor state of the constituent quark of color a and Nc = 3 
the number of quark colors. Since for a confining potential the strong interaction does not 
produce a pole at s = M 2 in the physical region (in the harmonic oscillator approximation of 
the quark model the Gaussian functions are smooth), the vertex function G v (s) can be related, 
as in Eq. ([8]), to a wave function representation of the form 

M") = Gv(s)/(s - M 2 ) = NpMs), (A7) 
where Np is a normalization factor and 



, / x 7T v s 2 — (rn% — m?) 2 1 . , . 

V2 yJS - (7712 - mi) 2 S d / 4 

In Eq. (|A8p . the function w(k) is chosen to be 

w(k) = exp (-4i/A; 2 /m 2 2 ) , (A9) 

where m\i is again the reduced mass. As in CLFD, we determine the normalization, Np, and 
fit the size parameter v so as to reproduce the experimental decay constants. In the dispersion 
approach the relativistic normalization Eq. @, by the appropriate choice for the wave function, 
reduces to the simple integral 

l = Nl w 2 (k)k 2 dk, with A^ P = — -( — =-) . (A10) 
Jo 7r 1 /4Vmf 2 / 

3. Decay constant of the pseudoscalar mesons 

According to the usual definition, the decay amplitude is 3' = (0| J 5/i |P) where J 5/i is the 
axial current. Since our formulation is explicitly covariant, we can decompose in terms of 
all momenta available in our system, i.e. the incoming meson momentum and u^. We have 
therefore: 

S^ = f P p fl + B^, (All) 

where fp is the physical decay constant. In an exact calculation of S„, B should be zero. Since 
uj 2 = 0, the decay constant can easily be obtained according to: 

fp = — • (A12) 
oj-p 

Using the diagrammatic rules of CLFD, we can calculate and including color factors, one 
gets, 



E^ = V3N P [ L>(x,0,R ± )Tr - 7 V(fe + m 2 )-^A {qq '\x, Ri) 7 5 (mi - fa) 

J{x,e,n ± ) L v2 



, (A13) 
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where the notation O is defined as usual by O = 7°0^7°. The decay constant is therefore given 
by: 



fp 



2V6N P \ 



D(x,6,H±) mi(l — x) + m 2 x 



A^ q Hx,-R 2 t 



(A14) 



Similarly, in the dispersion approach, taking into account soft rescatterings of constituent 
quarks, one obtains a series of dispersion graphs that involve the spectral density /9p(,s, mi, m 2 ) 
of the Feynman quark antiquark loop graph given in Eq. (|10p . These graphs yield the following 
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expression for the pseudoscalar decay constant 

fp = Np V Nc' dsm 1 + m2 



p P (s,m 1 ,m,2) (f)p(s). 



(A15) 



'(rai+m2) 2 77 S 

Then applying the normalization condition with the decay constants as constraints for 
modelling, we obtain the parameters listed in Table [DO 



APPENDIX B: DETAILS OF SPECTRAL DENSITIES 



Note that the double spectral densities for the pseudoscalar to scalar transition form factors 
A±(si, s 2 , g 2 ; mi, m 2 , m 3 ) in Eqs. (|72p and (|73p may be obtained from Melikhov 42| (section 
IIC) by the substitution mi into —mi. This substitution is the consequence of the different 
expressions of the operators in the trace entering in Eq. (j73[) in the case of p seudoscalar to 
scalar transition and in the case of pseudoscalar to pseudoscalar transition 
for completeness, we give here the explicit expression. One has 
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A±(si, s 2 , q 2 ; mi, m 2 , m 3 ) 



B±(si,s 2 ,q 2 
A(si,s 2 ,9 2 ) 



A(si,s 2 ,<? 2 ;mi,m2,m 3 ), 



(Bl) 



where 



B + (s 1 , s 2 , q 2 ) = b + (si, s 2 , q 2 ) [a(s 1 ,m 2 ,m 3 ) + a(s 2 , m 3 , -mi) - a(q 2 , -mi, m 2 )] 

+ a(q 2 , -m 1 ,m 2 )X(s 1 ,S2,q 2 ), (B2) 

B-(s 1 ,s 2 ,q 2 ) = b-(si,s 2 ,q 2 ) [a(s 1 ,m 2 ,m 3 ) + a(s 2 ,m 3 , -mi) - a(q 2 , -mi,m 2 )] 

+ [a(s 2 ,m 3 ,-mi) - a(si,m 2 ,m 3 )] A(si,s 2 ,g 2 ), (B3) 
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wi 



th 7 



6+(si, s 2 , q 2 ) = -q 2 {si + s 2 - q 2 + m\ + m| - 27713) ~ ( m i ~ m 2)( s i ~ s 2), 
6_(si,s 2 ,g 2 ) = (ml - 77i|)(2si + 2s 2 - q 2 ) 

+ {s\ - s 2 ){si + s 2 - q 2 + ml + ml - 2ml), 



(B4) 



(B5) 



(B6) 



with a(x, y,z) = x — (y — z) 2 . Furthermore, 

2 0(6^(si,s 2 ,9 2 ) - A(si,s 2 ,g 2 )A(g 2 ,m?,m!)) 

A(si,s 2 ,g ;mi,m 2 ,m 3 ) = TavU2~t T\ • 

16A 1 / 1 {s 1 ,s 2 ,q 2 ) 

The allowed intervals for the integration variables s% and s 2 are obtained by solving the step 
6>-function of Eq. (1B6L 



s 2 > (mi + m 3 ) 2 , 
«r( s 2 ; 9 2 ) < si < s^(s 2 ,g 2 ), 



(B7) 
(B8) 



with 



Si {s 2 ,q ) = 



q 2 ) + q 2 {m\ + to 2 ) 



(?77 2 — 777 2 )(?77 2 — m|) 



2mf 



A 1 / 2 (s 2 , 777 2 , m 2 )A 1 ' /2 ((7 2 , m 2 , 7773 
2m? 



(B9) 



The solution of the equation 



which reduces to 



, 9 + mf - m^ 2 2 

s 2 H == + "7 X - 777 3 = 0, 



(BIO) 



(Bll) 



so therefore the limit s 2 ((? ) appearing in Eq. (|74p is 



9 9 9 

g + 777J — 7772 

2^ 



+ 



7 2 + 777 2 — m| 



2\/Y 



+7773 — m|. 



(B12) 



Note that in Eqs. (|Blj) to (IB5h we have introduced, following Melikhov, a lightened writing for 
the functions B±(s\, s 2 , q 2 ) and b±(s\, s 2 , q 2 ) which, we stress, depend parametrically on the 
quark masses mi, m 2 and 1713. This is obviously the case also for sf(s 2 ,q 2 ) and s 2 (g 2 ). 



7 In the expression of 6_ (si, S2, q 2 ) given by Melikhov [4^1 (see his Eq. (2.76)) there is a misprint: the relative 
sign between the two term should be + as here in Eq. ()B5|I . 
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FIG. 1: Representation of the two-body wave function on the light front. 
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FIG. 2: The triangle diagram (leading contribution) and momentum flow in the weak-hadronic P — > 5 
transition amplitude of the CLFD approach. In the present case: ra-i > mi. 




FIG. 3: Same as in Fig. [2] but for the DR approach. 
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FIG. 4: The variation of x 2 /d.o.f. as a function of the mixing angle # m i X . It corresponds to the fit 1 
where 12 branching ratios are fitted with 2 parameters. The full and dashed lines correspond to the DR 
and CLFD results, respectively. 
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FIG. 5: Same as in Fig. H] but for the fit 2b where 9 branching ratios are fitted with 4 parameters. 
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FIG. 6: Effective form factors F (q 2 ) and Fi(q 2 ) [see Eq. (|75|)] calculated with the parameters of fit 2b 
for D — > /o(980) transitions. In the DR model, the full and dotted lines correspond to Fi(q 2 ) and Fo(q 2 ) 
respectively, and similarly for the dashed and dot-dashed lines in the CLFD model. 
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FIG. 7: Same as in Fig.[S]but for D s — > / (980) transitions [see Eq. (fBUj) for the definition of F (q 2 ) and 
Fi(q 2 )}- 
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FIG. 8: Same as in Fig. ©but for B ->■ /o(980) transitions. 
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FIG. 9: Same as in Fig. but for B s 



/o(980) transitions [see Eq. ([SO])]. 
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TABLE I: Fit 1 (12 branching ratios and 2 parameters): comparison of experimental with theoretical 
branching ratios. The fit parameters are found in Table ITT1 and a best fit yields x 2 /d.o.f. = 33.25/(12—2) = 
3.33 with DR and x 2 /d.o.f. = 30.63/(12 - 2) = 3.06 with CLFD. 







$mix 




N s 


CLFD 


(3.20 ±0.40) x 10^3 


32.0° ±4.8° 


(4.64 ±0.58) x 10^3 


2.00 


DR 


0.014 ±0.012 


41.3° ±5.5° 


0.021 ±0.017 


3.41 



TABLE II: The scalar-meson parameters, v n and 6> m i x , obtained in the CLFD and DR approaches with 
fit 1 (see Table QJ. Note that v s and Ns are given by Eq. (|75|) and by Eqs. (|2"6"|) or (|2"5j). respectively. 
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TABLE III: Fit 2a (12 branching ratios and 4 parameters): comparison of experimental with theoretical 
branching ratios. The fit parameters are found in Table ITVl and a best fit yields % 2 /d.o.f. = 12.73/(12 — 
4) = 1.59 with DR and x 2 /d.o.f. = 10.95/(12- 4) = 1.37 with CLFD. 
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v s 
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(3.19 ±0.32) x 10~ 3 
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0.23 ±0.05 


(4.62 ±0.46) x 10~ 3 


2.00 


DR 


0.014 ±0.012 
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0.35 ±0.42 


0.021 ±0.017 


3.42 



TABLE IV: Same as in Table U but annihilation included (fit 2a, see Table ITtT| . 
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TABLE V: Fit 2b (9 branching ratios and 4 parameters): comparison of experimental with theoretical 
branching ratios. The fit parameters are found in Table IVT1 and a best fit yields x 2 /d.o.f. = 6.82/(9 — 4) = 
1.36 with DR and x 2 /d.o.f. = 5.37/(9 - 4) = 1.07 with CLFD. 







^mix 
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v s 


N s 


CLFD 


(3.09 ±0.36) x 10~ 3 


31.5° ±5.0° 


0.23 ±0.25 


(4.49 ±0.52) x 10" 3 


1.93 


DR 


0.017 ±0.010 


41.6° ±7.1° 


0.34 ±0.41 


0.024 ±0.014 


3.84 



TABLE VI: Same as in Table HH but for fit 2b (see Table EJ. 
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q 2 


ml m 2 K m 2 p 




CLFD DR 
0.21 0.22 


CLFD DR 
0.28 0.18 


CLFD DR 
0.38 0.17 


F 1 D - tf0 (q 2 ) 


0.21 0.22 


0.33 0.22 


0.94 0.26 


F D '^ f0 (q 2 ) 


0.45 0.46 


0.67 0.41 


1.02 0.32 


F?^(q 2 ) 


0.45 0.46 


0.75 0.48 


1.86 0.53 



TABLE VII: Effective scalar and vector form factors F (q 2 ) and F\{q 2 ) [see Eqs. (|T9")) and ([50])] for various 
typical timelike momentum transfers, q 2 , in Dr s -\ — > /o(980) transitions in the CLFD and DR approaches, 
respectively (see Figs. [6] and [7]). 
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m 2 K 
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m 2 D 






CLFD 


DR 


CLFD 


DR 


CLFD 


DR 


CLFD 


DR 


F B ^ f "(q 2 ) 


0.12 


0.12 


0.13 


0.12 


0.14 


0.12 


0.23 


0.13 


F B->/ 0(g2) 


0.12 


0.12 


0.14 


0.12 


0.15 


0.13 


0.28 


0.15 




0.40 


0.29 


0.41 


0.30 


0.47 


0.30 


0.74 


0.29 


F^ h {q 2 ) 


0.40 


0.29 


0.43 


0.30 


0.51 


0.31 


0.89 


0.35 



TABLE VIII: Same as in Table EED but for B {s) -> / (980) transitions (see Figs. E and [9). 
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D, 


B 


B s 


CLFD 


N P 


9.976 


7.340 


5.880 


3.833 




V 


0.046 


0.061 


0.049 


0.059 


DR 


N P 


4.395 


3.443 


3.937 


2.646 




V 


0.043 


0.057 


0.049 


0.057 



TABLE IX: The pseudoscalar-meson parameters, Np and f, in the CLFD and DR approaches. The 
normalization Np is either calculated with Eq. (|A5[) or Eq. (|A10|) . The wave-function range parameter 
v, which enters the theoretical evaluation of the decay constants in Eq. (|A14[) and Eq. (|A15|) , is fitted to 
reproduce the experimental values of the decay constants of Eqs. ([50]) . 



